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SECTION  I 


INTRODUCTION 


Colloidal  suspensions  of  small  particles  in  electrolytes  are  known  to 
form  crystalline  arrays.  One  of  the  systems  receiving  the  greatest 
attention  is  that  of  monodisperse  latexes,  or  charged  polystyrene  spheres 
in  water  and  other  solvents  (References  1-3).  These  systems  arrange 
themselves  in  body-centered  cubic  (bcc)  (Figure  A-l)  or  face-centered  cubic 
(fee)  structures  (Figure  A-2)  depending  upon  the  particle  concentration 
(Reference  4).  These  crystalline  patterns  are  commonly  discerned  through 
ordinary  optical  microscopy,  as  illustrated  in  Figure  1.  There  is  also 
evidence  for  the  existence  of  hexagonal  close-packed  (hep)  structures 
(Figure  A-l)  under  certain  conditions  (Reference  5). 

These  arrays  are  sufficiently  well-ordered  to  produce  3ragg 
diffraction  of  light  in  the  optical  region  (References  5  and  6).  The 
measured  performance  of  a  very  high  optical  density  rejection  filter  made 
from  a  crysta 1  1  i ne  col  1 oida  1  array  of  polystyrene  spheres  has  recently  been 
reported  (Reference  7).  The  major  filter  data  from  this  reference  are 
summarized  in  the  following  table: 


TABLE  I 


_ Performance  Data  of  the  Filter  in  Reference  7 _ 

-  Minimum  Transmittance  <0.012 

-  Out  of  Band  Transmittance  >  902 

Bandwidth  (F.W.H.M.)  <  12  nm 

-  Wavelength  Shift  for  10°  Rotation  -4.4  nm:fc.0.3  nm 


*  (Obtained  from  Figure  2  of  Reference  7) 


We  have  successfully  explained  the  performance  data  of  this  new 
rejection  filter  using  dynamical  X-ray  diffraction  theory  (References  3  and  9) 
and  optical  scattering  theory  (Reference  10),  as  well  as  simple  considerations 
from  thin  film  optics  (Reference  11  and  12).  Extensive  expressions  are 
derived  for  the  bandwidth,  and  attenuation  of  the  filter,  in  which  case 
very  large  absorbance  values  are  predicted.  Some  of  our  work  in  shortened 
form  has  been  previously  presented  (References  13  and  14). 


SUCTION  11 


CRYSTAL  PARAMETERS 

The  crystal  line  col  loidal  array  filter  of  Flaugh  et  al. (Reference  7) 
was  made  with  a  rejection  wavelength  maximum  of  514.5  nm.  We  shall 
designate  this  wavelength  maximum  in  air  as  J\sl,  which  is  related  to  the 
di stance  d  between  the  planes  of  the  crystalline  array  by  the  Bragg 
condi tion 

ii  A*  =  2nd  si  n  8  .  ( 1 ) 

Here,  M  is  the  diffraction  order,  n  is  the  refractive  index  of  the 
colloidal  suspension,  and  8  is  the  usual  Bragg  angle  within  the  suspension. 
The  geometric  and  parametric  relationships  for  the  incident  and  diffracted 
light  rays  are  depicted  in  Figure  2.  Taking  the  refractive  index  to  be 
equal  to  1.333,  that  of  water,  we  obtain  d  =  193.0  nm,  for  first-order 
diffraction.  Flaugh  et  al.  state  that  the  filter  crystal  structure  was 
face-centered  cubic,  and  the  filter  planes  are  therefore  (111)  planes.  The 
fee  unit  cell  is  i  1  lustra  ted  in  Figure  3,  from  which  we  may  evaluate  the 
lattice  parameter  as  a  =  dVT,  or  334.3  nm.  Our  values  for  d  and  a_  do  not 
agree  with  those  of  Flaugh  et  al.  because  they  neglected  the  refractive 
index  of  the  liquid,  and  in  the  case  of  a,  they  seem  to  have  a 
typograph ica 1  error. 

The  diameter  JQ  of  the  polystyrene  spheres  used  by  Flaugh  et  al. 
was  91  nm.  Ihe  volume  fraction  occupied  by  the  spheres  is  therefore 
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or  4.22  percent. 

As  found  experimentally  by  lliltner  et  al.  (References  2  and  6),  the 
refractive  index  varies  according  to 

n  =  1.333  +  0.27  0  (; 


with  the  result  that  n  =  1.344,  for  the  present  case.  Thus,  there  is  only 
a  slight  error  in  using  the  refractive  index  of  water  in  calculating  d  and  a 

Alternatively,  one  may  combine  Eqs.(l-3)  and  solve  the  resulting  cubic 
equation  for  d,  yielding  d  =  191.3  nm.  This  is  close  enough  to  the  result 
of  the  simplified  calculation  for  d,  that  our  first  result  is  sufficient 
for  use  in  the  rest  of  our  calculations. 

Another  parameter  we  shall  need  is  D,  the  next  nearest  neighbor 
spacing  (Figure  4).  This  is  obtained  by  D  =  a  /  ~V2~,  or  236.4  nm.  Also,  the 
colloidal  array  thickness  was  3.97  x  10^  nm,  equivalent  to  2057  (111) 
planes  (Figure  5). 


Figure  5.  Macroscopic  Arrangement  of  the  Colloidal 
Array  Filter.  A^  and  A„  are  Plexiglass 
Plates;  C  is  the  Colloidal  Array  Held  in 
Place  by  a  Teflon  Spacer. 


SECTION  III 


ANGLE  SHIFT  OF  THE  REJECTION  BAND  MAXIMUM 


The  angle  shift  of  the  reflection  band  maximum  is  shown  in  Figure  6. 
If  we  designate  S A&  to  be  the  wavelength  shift  of  the  rejection  band 
maximum  upon  rotation  of  the  filter  by  an  angl e ^rel ati ve  to  the  incident 
radiation  (See  also  Figure  5),  application  of  Snell's  Law  produces  the 
following  relation  between  &Ao  and  known  quantities: 


f  = 


For  small  values  of 


- 


this  expression  reduces  to 


(4) 


A0 


ia2 

znz 


(5) 


a  result  well-known  in  interference  filter  theory  (References  11  and  12). 

f  laugh  et  al.  state  that  Q0i  s  506.7  nm  minus  514.5  nm,  or  -7.8  nin,  for 
an  angle  of  incidence  of  10°.  However,  careful  measurement  of  the  data  from 
their  Fig.  2  yields  a  value  of  &Ao  =  -4.4  i  0.3  nm,  where  the  error 
indicated  is  the  precision  in  reading  SAo  from  their  figure.  If  we  use 
n  -  1.344  ,  we  calculate  &A0  to  be  -4.3  nm  using  Eq.  4,  consistent  with  the 
value  we  measure  from  their  data.  We  believe  the  stated  wavelength  shift 
of  I  laugh  et  al.  to  be  in  error  because  it  deviates  greatly  from  the 


theoretical  value. 


SUCTION  IV 


BANDWIDTH 


By  bandwidth,  we  shall  use  the  definition  of  the  full  width  of  the 
linear  transmission  band  at  one-half  the  maximum  intensity  (F.W.H.M.) 
which  we  shal  1  1 abel  aA0.  To  calculate  the  bandwidth,  we  must  employ  the 
dynamical  theory  of  X-ray  diffraction  (References  8  and  9).  The  necessity 
for  this  approach  is  readily  seen  by  referring  to  Figure  7.  The  incident 
beam  of  intensity  IQ  loses  a  small  fraction  of  its  energy  by  reflection 
from  each  successive  lattice  plane.  In  addition,  the  initially  reflected 
beam  is  reflected  a  second  time  to  a  direction  parallel  to  the  incident 
beam.  These  effects  are  neglected  in  the  kinematic  theory  of  x-ray 
diffraction  usually  used  by  crysta  1 1  ographers  to  determine  the  structure  of 
ina  tter . 

We  shall  make  use  of  the  development  of  this  theory  by  Ewald  as  later 
expounded  by  Zachariasen  (Reference  8).  The  starting  point  shall  be  his 
Eq.  (3.93)  relating  the  dielectric  constant  €  (?)  and  the  polarizability  per 
unit  volume  (r*),  both  a  function  of  the  lattice  position  vector  7f: 

e(r)  =  i-v4Tr^^)  .  (6) 


VJe  nay  further  express  (r)  as  the  product  of  <Xs  ,  the  pol  ari  zabi  1  i  ty  of 

a  single  sphere,  and  N(r),  the  sphere  spatial  distribution.  The  sphere 
pol ari zabi 1 i ty  is: 


s 


4*7 f(V*l+2.) 


Vs 


(/) 


where  m  is  the  ratio  of  the  refractive  index  of  the  spheres  to  that  of  the 
surrounding  liquid  medium,  and  Vs  is  the  volume  of  a  sphere  (Reference  10). 
The  quantity  4  ? Tot  (?)  may  be  replaced  by  a  new  function^  (r),  which,  being 
periodic,  can  be  expanded  as  a  Fourier  series: 

'+V)  =  4ir=<('f?)  =  5:g;  .  (3) 

H 

In  this  equation!^  is  the  reciprocal  lattice  vector  satisfying  the  Laue 
equation 


(9) 


and  k,  and  k2  are  the  incident  and  scattered  wave  vectors.  It  can  then  be 
shown  that 


%  =  4-?ft£k..  ,  no) 

where  F^  is  the  crystal  structure  factor  when  the  Laue  equation  is 

3 

sa  tisfied,  and  V  is  the  vol ume  of  a  crystalline  unit  cell,  equa 1  to  a 
for  the  fee  lattice. 

For  waves  diffracted  off  (111)  planes  in  the  fee  lattice,  FH  is 
evaluated  to  be  4f,  f  being  the  scattering  factor  for  a  sphere.  To 
evaluate  _f,  we  make  use  of  the  scattering  theory  of  small  spheres  developed 
by  Van  De  Hulst  (Reference  10).  With  the  definition 

,  (in 


14 


the  condition  for  ordinary  Rayleigh  scattering  is  (Reference  10) 


X«  1  ,  (12) 

which  is  true  of  Region  61  of  Figure  3  and  Table  II.  If  Eq.  (11)  does  not 
hold,  some  researchers  believe  that  Mie  scattering  theory  must  be 
immediately  utilized.  However,  in  Region  1,  one  is  yet  able  to  use  a 
modified  Rayleigh  theory,  known  as  Rayl ei gh-Gans  theory,  provided 

x(Vn-l)«l  .  (13) 

To  further  understand  the  application  of  the  above  conditions  to  the  case 
of  the  colloidal  array  filter,  we  have  plotted  x(m-l)  and  x  vs.  Do/d 
in  Figure  9.  For  our  case.  Do/d  is  0.472  and  x  is  0.73,  which  would 
be  too  large  for  ordinary  Rayleigh  scattering  theory.  Additional  justification 
for  using  Rayl eigh-Gans  theory  may  be  found  in  Reference  10,  Sections  18.4 
and  19.12. 

Continuing  to  follow  the  development  of  Van  De  Hu  1  s t ,  we  obtain 

P  =  G(i\)  -  ~r("sin  ^  -  a  cos  u.)  j  (14) 

where 

i  i  _  2JL  c  |'in  0 

u  -  ^  b  m  o  *  (15) 


L5 


V  .  V  V 


Figure  8.  Regions  of  Validity  of  Various  Scattering 
Theories  in  the  m-x  Domain  (After 
Reference  10). 
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In  this  equation,  _A_  is  the  wavelength  of  light  in  the  colloidal 
suspension,  equal  to  A„/n.  Note  that  0,  the  Bragg  angle,  is  one-half 
the  "scattering  angle"  usually  defined  in  optical  scattering  (Reference  10) 
Combining  Eqs.  1,  7,  10,  14,  and  15,  we  find 

(s'n  a  'acoSU)  no, 

where 


(17) 


Zachariesen  (Reference  8)  has  shown  that  for  a  thick  crystal  and  no 
absorption,  the  bandwidth  is 

wY  1  K  %  [  t 

Yb”  sin*  6 


(18) 


where  wy  is  a  numerical  factor  2/3  VI  =  1.155  for  Ewald  theory  and  b  is 
t!ie  ratio  of  direction  cosines  of  the  light  rays,  which  is  unity  for  Bragg 
diffraction.  The  polarization  factor  K  is  unity  for  polarization  and 

'.os  2  9 1  f  o  r  7f  po  1  ar  i  za  t  i  on.  The  correctness  of  using  X,  instead  of  A 
in  this  equation  is  readily  a  seer ta i nab  1 e  from  Eq.  1,  and  the  expression 
of  Cuchariesen  for  A_e,  the  angular  bandwidth.  Inserting  Eq.  16  into 
l q .  18  yields 


(sin  a  -  u,  cos  u) 

Sm2-  © 


.  (19) 


fins  expression  resembles  the  bandwidth  A  8  on  the  glancing  angle 
scale  derived  by  other  authors  (References  2,  6,  and  15),  but  with  some 

Ihe  /\0/2  and  sin^Q  terms  are  unique  to  (  aAj',  and 


important  differences. 


the  previous  derivations  have  not  considered  polarization.  Hiltner  et  al. 
'References  2  and  6),  used  a  value  of  w  =  3/4  1.061 ,  fo  1  1 owi ng  the 

L/arwin  theory  of  dynamical  diffraction  (Reference  8).  They  also  used  the 
approx  i,na  t  ion  in  the  pol  arizibi  1  i  ty  that  (m^  +  2)^3.  .‘-lei  one  and 
Rustichelli  (Reference  15)  used  the  approximations  _f  =  1  and  =  1. 

Eq.  19  is  accurate  for  the  bandwidth  if  the  col  1  oidal -ai r  interface 
is  neglected.  A  correction  based  on  Snell's  law  must  be  included  for  the 
bandwidth  measured  in  air.  If  nQ  is  the  refractive  index  of  air  and  8Q 
is  the  glancing  angle  in  air,  it  can  be  shown  that  the  corrected  bandwidth 
i  s 

aX0 

where 


9  Tr2- 


si  n0 


Jsing  m  =  1.60/1.333  and  other  values  previously  stated,  we  calculate 
aAp  =  10.3  nm.  This  is  in  excellent  agreement  with  the<12  nm 
bandwidth  claimed  by  Flaugh  et  al.  (Reference  7).  Reasons  for  lack  of 
closer  agreement  between  theory  and  the  measured  value  may  include 
absorption,  crystalline  mosaicity,  lattice  distortions,  and  stacking 
fau Its. 

It  is  appropriate  to  consider  how  to  design  a  filter  with  a 


tailorable  bandwidth,  from  Eq.  20,  we  can  see  that  there  is  some  freedom 
in  choosing  the  refractive  index  of  the  spherical  particles  or  of  the  host 
medium.  The  number  of  materials  for  the  spheres  is  quite  limited,  but 


rrr.  yrr 
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there  are  several  liquids  known  to  support  crystalline  structures 
(Reference  2).  However,  the  refractive  indices  of  all  these  liquids  are 
greater  than  that  of  water,  which  would  only  broaden  a  filter  rejection 
band.  We  are  thus  led  to  consider  the  effect  of  smaller  particle 
diameters  upon  the  filter  bandwidth.  We  shall  define 

H(u)  3  (sinu-ucosu)  (22) 

and  plot  this  as  a  function  of  Do/d  in  Figure  10.  This  graph  demonstrates 
that  a  narrower  bandwidth  is  most  easily  obtained  by  using  smaller 
diameter  polystyrene  spheres. 

For  the  crystal  1 ine  colloidal  array  filter  of  Flaugh  et  a  1 . 

(Reference  7),  u  was  1.48,  corresponding  to  H(u)  =  0.864.  In  this 
region,  it  can  be  readily  seen  from  Figure  10  that  the  cubic  approximation 
for  H(u)  should  not  be  used. 


SECTION  V 


ATTENUATION 


Zachariasen  (Reference  8)  has  shown  that  the  maximum  reflectance  R^( 
for  3ragg  diffraction  takes  the  form 

R^o  -  (tanh  A)2  .  (i 

In  the  absence  of  absorption, 


y  k  1  Vh  (  t. 


A  X 


(; 


where  t0  is  the  depth  of  the  wave  into  the  crystal  along  a  normal  to  the 
crystal  surface,  and  To  =  sin  6.  Combining  Eqs.  16  and  24  gives 


.  _  2. 
ir( 


^-acosa)  17VT 


.  (; 


ima+2. JV  J  sind 

The  transmittance  TA#  at  the  wavelength  Aj,  may  then  be  expressed  as 


t;  =  l- R.- 


cos 


Hf) 


1  -l 


(; 


where 


_  WyfAo)2  sin  $ 

z  ~  V  fa  Ac) 


For  normal  incidence,  t  is  calculated  to  be  9.42  x  ICr  nm  for  the 
colloidal  filter  of  Flaugh  et  al.  This  value  of  t  is  equivalent  to  48.8 
layers  which  may  be  compared  to  the  total  thickness  of  Flaugh's  filter 
which  is  2057  layers.  The  transmittance  may  further  be  configured  as 

Tx.=  cosli(r) 

where  L  =  48.8  and  2  is  the  numbers  of  layers  penetrated  by  the  incident 
beam.  Eq.  28  is  plotted  in  Fig.  11,  where  the  results  may  be  seen  to  be 
consistent  with  the  experiments  of  Flaugh  et  al.  who  could  measure 
transmittance  values  only  as  low  as  10-4.  Also,  the  theory  demonstrates 
that  extremely  large  absorbance  values  are  possible  in  colloidal  filters 
much  thinner  than  that  of  Flaugh  et  al . 

Since  the  time  of  our  original  reports  (References  13  and  14),  Prof. 
Asher  (Reference  15)  has  measured  optical  densities  greater  than  eleven, 
in  the  direction  of  extremely  large  attenuation  predicted  by  our  theory. 
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SECTION  VI 


SUMMAR V 


The  dynamical  theory  of  X-ray  diffraction  and  light  scattering  theory 
have  been  applied  to  the  calculation  of  the  bandwidth  and  attenuation  of  a 
crystal  1  ine  col  1 oida  1  array  optical  filter.  The  1 imi ts  of  the  scatteri ng 
theory  dependent  upon  particle  diameter  and  refractive  index  have  been 
thoroughly  defined.  The  theoretical  bandwidth  is  in  excellent  agreement 
with  the  previously  published  experimental  value.  Design  parameters  for 
narrower  bandwidth  filters  are  discussed. 

The  theoretical  attenuation  results  are  consistent  with  measured 
transmittance  values.  The  theoretical  attenuation  function  predicts  that 
actual  filters  may  be  constructed  that  are  thinner  and  more  efficient  than 
the  first  experimental  ones,  as  has  been  substantiated  by  the  most  recent 
transmi ttance  measuranents. 

The  theoretical  angle  shift  of  the  rejection  band  maximum  has  also 
been  calculated,  and  is  in  excellent  agreement  with  experimental  data. 
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APPENDIX  A 


ADDITIONAL  CRYSTAL  STRUCTURES 


hcc 

Hexagonal  close-packed 


fee 

Faced-centered  cubic 


Figure'  A—  1  .  Cl ose-Packod  Crystal  Structures.  The 

(111)  Axis  is  in  the  Vertical  Direction. 


